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A  Two-Stream,  Multilayer, 

Spectral  Radiative  Transfer  Model  for  Sea  Ice 


DONALD  K.  PEROVICH 


INTRODUCTION 

The  reflection,  absorption,  and  transmission  of  shortwave  radiation  by  snow  and  ice 
covers  are  significant  in  a  variety  of  geophysical  problems.  Spectral  reflectance  is  important 
in  the  interpretation  of  imagery  from  remote  sensing  instruments  operating  at  visible  and 
near-infrared  wavelengths.  Maykut  and  Untersteiner  (1971),  Grenfell  and  Maykut  (1977), 
and  Maykut  and  Perovich  (1987)  have  established  the  critical  role  played  bv  shortwave  radi¬ 
ation  in  the  heat  balance  of  a  sea  ice  cover.  Primary  productivity  and  other  biological  activity 
under  a  sea  ice  cover  are  strongly  affected  by  the  amount  and  spectral  composition  of  trans¬ 
mitted  light  (SooHoo  et  al.  1987). 

Because  of  this  established  importance,  there  is  a  substantial  data  base  of  optical  meas¬ 
urements  of  reflected  and  transmitted  light  in  snow  and  ice  covers.  Observational  data  on  the 
optical  properties  of  snow  and  of  sea  ice  are  reviewed  in  Warren  (1982)  and  Perovich  et  al. 
(1986),  respectively.  These  observations  indicate  that  understanding  radiative  transfer  in 
snow  and  ice  is  considerably  complicated  by  the  large  temporal  and  spatial  variability  in  the 
physical,  and  therefore  optical,  properties  of  the  medium.  For  example,  observations  indi¬ 
cate  that  as  a  snow  cover  ages,  and  metamorphism  increases  its  grain  size,  there  is  a  decrease 
in  albedo  (Grenfell  and  Perovich  1981).  The  physical  structure  of  sea  ice  can  be  even  more 
complex,  with  vertical  variations  in  vapor  inclusions,  ice  temperature,  salinity,  and  brine 
volume.  Observations  have  demonstrated  that  the  optical  properties  of  sea  ice  depend  on  the 
ice  temperature,  brine  volume  (Perovich  and  Grenfell  1981 )  and  surface  conditions  (Grenfell 
and  Maykut  1977). 

This  complexity  necessitates  use  of  theoretical  models  both  to  interpret  and  to  extrapolate 
the  observational  data.  Over  the  years  several  different  models  have  been  formulated  to  de¬ 
scribe  radiative  transfer  in  snow  and  ice.  Warren  (1982)  provides  an  excellent  summary  of 
radiative  transfer  models  that  have  been  applied  to  snow.  One,  in  particular,  developed  by 
Wiscombe  and  Warren  (1980),  has  proved  very  useful  for  snow  applications.  Mullen  and 
Warren  (1988)  used  a  modified  delta-Eddington  approach  to  model  spectral  albedos  of  lake 
ice.  Perovich  and  Grenfell  (1982)  and  Grenfell  (1983)  adapted  the  discrete  ordinates  method 
of  Chandrasekhar  (1960)  to  sea  ice.  The  Dunkle  and  Bevans  (1956)  model,  based  on  Schuster 
(1905),  calculates  the  upwellingand  downwelling  irradiance  in  a  snowpack.  The  model  was 
used  by  Grentell  and  Maykut  (1977),  Grenfell  (1979),  and  Perovich  and  Grenfell  (1981)  to 
interpret  field  and  laboratory  measurements  of  the  optical  properties  of  sea  ice. 

These  models  all  have  features  that  both  recommend  and  limit  them  for  particular  appli¬ 
cations.  Often  there  are  tradeoffs  between  detailed  physics  and  computational  simplicity.  In 
this  report  we  will  develop  and  present  a  two-stream,  multilayer  model  of  radiative  transfer 
at  visible  and  near-infrared  wavelengths.  While  this  discussion  is  directed  specifically  to¬ 
ward  a  sea  ice  cover,  the  formulation  is  sufficiently  general  to  be  easily  adapted  to  other  tran- 
lucent  media  such  as  lake  ice,  river  ice,  or  turbid  water.  The  major  limitation  of  this  model 


lies  in  its  casual  treatment  of  the  physics  of  scattering,  in  particular  in  its  assumption  of 
isotropic  scattering  and  the  lack  of  information  concerning  the  angular  distiibution  of  the 
radiation  field.  However,  since  our  interest  is  in  irradiance,  the  radiance  integrated  over 
angle,  this  is  not  a  major  difficulty.  The  model  has  three  primary-  advantages.  Computationally 
it  is  simple  and  fast,  avoiding  the  numerical  difficulties  of  more  advanced  methods.  Thus  it 
can  be  run  quickly  on  a  personal  computer,  allowing  a  wide  range  of  cases  to  be  investigated 
easily.  In  addition,  this  two-stream  model  can  directly  exploit  the  large  observational  data 
base  of  sea  ice  optical  properties  compiled  bv  Grenfell  and  Mavkut  (1977),  Grenfell  (1979), 
and  Perovich  and  Grenfell  (1981).  Finally,  only  a  qualitative  depiction  of  the  ice  structure 
(blue  ice,  white  ice,  melting  ice)  is  needed,  rather  than  the  detailed  statistical  description  of 
the  scattering  inhomogeneities  demanded  by  more  sophisticated  models. 


THEORY 

The  model  is  of  two  streams  in  that  it  calculates  a  downvvelling  (F;)  and  an  upwelling  (F.) 
irradiance  in  the  medium.  The  irradiance  is  defined  as  the  total  energy  per  unit  area  per  unit 
time  ( V V /  nv)  incident  on  a  horizontal  plane  from  the  lower  (F.)  and  upper  (F  )  hemispheres. 
In  the  medium  each  layer  is  considered  to  be  plane  parallel,  i.e.  homogeneous  in  z  and  infinite 
in  extent  in  .v  and  y.  The  optical  properties  of  a  layer  are  defined  in  terms  of  wavelength-de¬ 
pendent  scattering  (r )  and  extinction  (k. )  coefficients.  The  downvvelling  and  upwelling  ir- 
radiances  are  governed  by  the  coupled  first-order  differential  equations: 


dF  (z,X)  = 

-k.fFdz,X)dz 

-r.F(z.X)dz 

+r.Fr(z,X)dz 

Change  in 
downwelling 

Loss  to 
absorption 

Scattering  loss 
from  downvvelling 

Scattering  gain 
from  upwelling 

dF.(z,X)  = 

kF.(z,X)dz 

-rF,(  z,X)dz 

+r.F^z,X)dz 

Change  in 
upwelling 

Loss  to 
absorption 

Scattering  gain 
from  downwelling 

Scattering  loss 
from  upwelling 

where  r  is  the  depth  within  the  medium  (increasing  downward)  and  X  is  the  wavelength  of 
the  light.  Scattering  (r )  and  absorption  (k. )  are  treated  through  wavelength-dependent  coef¬ 
ficients.  Together  the  scattering  and  absorption  coefficients  define  the  extinction  coefficient 
k.  where  k.  =  V(A‘. :  +  2k.  rt).  If  there  is  no  scattering  eq  1  reduces  to 

dFt(z,X)  =  -k.Fiz,X)dz, 

which  when  integrated  yields  the  familiar  Bouguer-Lambert  law 
Ft(:,)v)  «  e~*'-'  . 

Equation  2  reduces  to  F  .(z,A)  =  0  for  the  no-scattering  case.  Equations  1  and  2  can  be  reformu¬ 
lated  as  uncoupled  second-order  differential  equations  (see  App.  A  for  details)  with  general 
solutions  of 


F,(z,X)  =  A  sinh(tcz)  +  B  cosh(icz)  (4) 

F.(z,X)  =  C  sinh(xzz)  +  D  cosh(icz).  (5) 


2 


This  method  converges  to  a  simple  expo¬ 
nential  d^cav  tor  optically  thick  media,  but 
as  Grenfell  (1*579)  points  out  it  offers  sub¬ 
stantial  improvements  in  accuracy  over 
exponential  decav  for  opticallv  thin  cases. 

The  geometry  of  the  multilayer  model  is 
schematically  illustrated  in  Figure  1.  We 
consider  N  layers,  with  each  layer  i  having 
thickness  H  and  optical  properties  defined 
bv  the  constants  Ay  and  r.  .  YVe  assume  that 
there  is  no  specular  reflection  at  the  inter¬ 
faces  between  layers.  To  specify  the  solution 
to  eq  4  and  5  for  a  layer  it  is  necessary  to  de¬ 
termine  the  four  constants  .4,  B,  C,  and  D. 
Hence  42V  equations  are  needed  to  solve  an 
.Y-laver  system. 

Two  equations  are  provided  by  the  surface  a 
face  we  have 


Incident  irradiance  (F0 ) 


(l  -  R 0  )  F„  (0,  X )  I  4  Specular  Reflection  (R0F0 ) 


t 


O  'No  upwetlmg  radiation  from  the  water  column) 

;hjv  2.  Schematic  of  tzco-strenm,  n-lai/cr  model. 
id  bottom  boundary  conditions.  At  the  sur- 


F  ,UU)  =  (l-R,i)F.+RiiF; ,(0,X),  (6) 

where  Fn  is  the  incident  solar  irradiance  and  R ,  is  the  specular  reflection  at  the  surface.  At  the 
bottom  boundary  we  assume  that  there  is  no  specular  reflection  and  there  is  no  upwelling 
from  below,  yielding  a  boundary  condition  of 

F-n(Hv/.)  =  0.  (7) 

YVe  require  continuity  at  the  interfaces  between  layers.  These  matching  conditions  provide 
2(N-1)  equations  of  the  form 

F  jH,,A)  =F  ,..(0,7.)  (8) 

F.(H,,W  =FIiM(0,A).  (9) 

In  eq  6  through  9  the  two  subscripts  on  F  define  the  direction  (i,T)  and  the  layer  number  i, 
and  H  is  the  layer  thickness.  The  remaining  2 N  equations  are  obtained  by  applying  the  di¬ 
fferential  equations  1  and  2  at  the  top  of  every  layer.  This  is  done  at  the  top  of  the  layer,  since 
c  =  0  simplifies  the  resulting  equations. 

Substituting  the  general  solutions  (eq4and  5)  into  these  4iV  equations  yields  the  following 
system  of  equations  for  each  wavelength. 


(1  -  R  )F  +  R  D  =  B. 

«>  o  * »  I  I 

CN  sinh(KNHN)  +DN  cosh(K'NHN)  -  0 
A  sinh(K  H)  +B  cosh (kH)  =  B  , 

C  sinh(K  H  )  +D  cosh(K  H )  =  D  . 

;  I  i  i  it  I*-! 

Alt  +  B(t+  r)-D  r  =  0 

II  il  I  II 

dr  +  Ck  -  (Ai  +  r  )D  =  0 


(top  boundary  condition) 

(bottom  boundary  condition) 

(iV-1  matching  conditions,  2=1,  N- 1) 
(N-l  matching  conditions,  i  =  1,  N-l) 
(N  equations  for  i  =  1  ,N) 

(N  equations  for  i  =  1 ,2V) . 
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Figure  2.  Matrix  representing  si/stein  of  equations  for  tuv-streain,  n-lni/er  model.  chdesiguntes  a  hyperbolic  cosine  nmi 
sh  denote s  a  In/perbolic  fine. 


Figure  2  presents  the  system  of  equations  for  an  N-laver  medium  in  matrix  form.  Once  the 
system  is  solved  and  values  of  A ,  B,  C  and  D  determined,  the  upwelling  and  downwelling 
irradiance  can  be  calculated  anywhere  in  the  medium.  When  computing  irradiances  the 
depth  ;  is  with  respect  to  the  layer  rather  than  the  entire  medium.  This  coi.ventii  n  improves 
the  numerical  accuracy  of  the  solution. 

Two  parameters  of  particular  interest  are  the  spectral  albedo  (cr,  and  the  transmittance 
CT  ).  The  spectral  albedo  is  the  ratio  of  the  reflected  and  incident  irradiance,  or 

a  =  (R/.+D,)  /  F„ 

The  spectral  transmittance  is  the  ratio  of  the  transmitted  and  the  incident  irradiance,  or 
T.  =  [AN  sinh(x\HN)  +  cosh(tcN.HN)]/Fi. 

These  spectral  values  can  be  integrated  over  wavelength  to  give  bulk  values  of  albedo  (au) 
a„  =  jct  F,, (X)d/ /\F{)(X)  dX 
and  transmittance  (Tu) 

T0  =  \T,F[tO i)  dX/iFu{X)  dX. 

For  biological  purposes,  the  downwelling  spectral  irradiance  at  the  underside  of  the  ice  is 
of  considerable  interest.  Since  it  is  the  number  of  photons  of  particular  wavelengths  that 
controls  biological  activity,  irradiance  is  presented  in  terms  of  photons  of  wavelength  X  per 
unit  area  per  time,  rather  than  the  more  familiar  units  of  W/irr  (SooHoo  et  al.  1987).  The  basic 
unit  for  this  is  an  einstein  (E),  which  is  equal  to  the  energy  of  a  mole  of  photons  of  wavelength 
More  precisely  an  einstein  is  defined  as 

E  =  nhc/X, 

where  n  =  Avogadro's  number  (6.025  x  1021),  h  =  Planck's  constant  (6.625  x  10_2J  J  s),  and  c 
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=  speed  of  light  (3  x  lO1'  run/ s).  The  total 
i  photosvnthetically  active  radiation  (PAR) 

7  is  equal  to  the  irradiance  (units  of  mE/ m: 
s)  integrated  from  400  to  700  nm. 

r* 

S 

THE  SEA  ICE  MODEL 

To  illustrate  the  efficacv  of  the  multi¬ 
layer  formulation,  we  now  apply  it  to  the 
specific  case  of  a  sea  ice  cover.  The  input 
parameters  needed  for  the  model  are  spec¬ 
tral  scattering  and  extinction  coefficients 
for  each  medium  type  considered,  layer 
thickness,  and  optionally  the  incident 
spectral  irradiance. 

For  the  sea  ice  model  nine  medium 
types  were  included:  1)  cold  dry  snow,  2) 
4  melting  wet  snow,  3)  ice  colder  than  the 
eutectic  point,  4)  surface  scattering  layer 
d  of  white  ice,  5)  the  interior  portion  of 

3  white  ice,  6)  cold  blue  ice,  7)  melting  blue 

ice,  8)  bubble-free  fresh  ice,  and  9)  clear 


waveiengf 


Figure  3.  Sped  nil  extinction  coefficients  for  nine 
medium  topes  used  in  tlu  sen  ice  radiative  transfer 
model;  1  >  dry  snow,  2)  melting  snow,  3)  ice  colder  (1979)  (Smith  and  Baker  1981).  The  melting 
than  the  eutectic  point,  -!■)  surface  scattering  layer  of  wet  snow  and  cold  blue  ice  values  were 
white  ice,  5)  interior  portion  of  icliite  ice,  6)  cold  blue  obtained  by  adjusting  the  scattering  coef- 
ice,  7)  melting  blue  ice,  S)  bubble-  free  fresh  ice,  and  ficients  for  cold  drv  snow  and  melting 
9)  clear  arctic  water.  blue  ice  in  order  to  reproduce  spectral  al¬ 

bedos  reported  in  Grenfell  and  Perovich 
(1984)  for  melting  snow  and  cold  blue  ice.  Extinction  coefficients  beyond  800  nm  have  been 
extrapolated  from  observations  and  are  somewhat  suspect.  Because  of  experimental  errors 
due  to  light  leakage  and  lack  of  spectral  resolution,  extinction 
coefficients  at  longer  wavelengths  are  quite  possibly  larger 

than  reported  here.  Table  L  Scatteri"g  coeffi" 

Values  for  the  scattering  coefficients  were  not  explicitly  c'ents  )  f°r  n*ne  medi- 
reported  in  any  of  the  above  papers,  though  spectral  albedos  um  types. 

wrere  presented.  Therefore  we  determined  scattering  coeffi-  Medium  'coefficient 

cients  by  entering  the  extinction  coefficients  from  Figure  3 

into  the  two-stream  model  and  varying  the  scattering  coef-  Cold  dry  snow  800 

ficient  to  reproduce  the  reported  albedos.  In  the  visible  and  Melting  wet  snow  lbii 

.  .  ,  .  ...  .  .  .  .  Ice  buow  eutectic  ml) 

near  infrared,  the  scattering  inhomogeneities  are  much  larg-  Vvhite  ice  scattering  120 

er  than  the  wavelength,  and  the  scattering  coefficient  does  white  ice  interior  2.5 

not  vary  with  wavelength  (Grenfell  1983,  Bohr^n  and  Huff-  Cold  blue  ice  1.8 

man  1983).  The  bubble-free  ice  and  clear  arctic  water  cases  1L°  1 ~ 

were  assumed  to  be  purely  absorbing  media,  Scattering  co-  C)car  arctjc  watcr  ()  0 
efficients  (Table  1 )  varied  from  zero  for  bubble-free  fresh  ice  Note:  Xhc  scdttorillg  aicffidonts 
and  arctic  water  to  800  nr1  for  cold  drv  snow  arc  constant  with  wavelength. 


4  arctic  water.  Spectral  extinction  coef- 
i  ficients  for  these  nine  types  are  plotted  in 
toco  Figure  3.  The  extinction  coefficients  for  ice 
types  1,  4,  3,  and  7  came  directly  from 
Grenfell  and  Mavkut  (1977),  type  3  from 
me  Perovich  (1979),  type  8  from  Grenfell  and 
ter  Perovich  (1  381 )  and  type  9  from  Grenfell 


Table  1.  Scattering  coeffi¬ 
cients  (nr1)  for  nine  medi¬ 
um  types. 


Suiflenng 
Medium  inefficient 

Cold  drv  snow  800 

Melting  wet  snow  160 

Ice  be  iow  eutectic  IhO 

White  ice  scattering  120 

White  ice  interior  2.5 

Cold  blue  ice  1.8 

Melting  blue  ice  1.2 

Bubble-free  ice  0.0 

Clear  arctic  water  0.0 

Note:  The  scattering  coefficients 
arc  constant  with  wavelength. 
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Figure  4.  Normalized  values  of  spectral  incident  irradiance  under  arctic  con¬ 
ditions  of  1 )  clear  skies,  2)  light  clouds,  solar  disk  visible,  3)  cloudy,  solar  disk 
barely  visible,  and  4)  Heavy  overcast,  solar  disk  not  visible  ( from  Grenfell  and 
Perovichl984).  These  values  are  multiplied  by  wavelength  integrated  incident 
irradiance  (400-2400  nm)  to  obtain  spectral  incident  irradiances. 


Spectral  albedos  and  transmittances  do  not  depend  on  the  spectral  composition  of  the 
incident  radiation,  so  if  these  are  the  parameters  of  interest  FJfk)  is  not  needed.  However, 
there  are  other  situations,  particularly  in  biological  applications,  where  absolute  values  of 
upwelling  or  downwelling  irradiance  are  desired,  in  these  cases  FU(X)  must  be  specified. 
Representative  arctic  spectra  of  incident  solar  irradiances  measured  at  Point  Barrow,  Alaska, 
were  incorporated  into  the  model  (Grenfell  and  Perovich  1984).  The  sky  conditions  for  the 
four  curves  presented  in  Figure  4  were  1 )  clear  skies,  2)  light  clouds,  solar  disk  clearly  visible, 
3)  cloudy,  solar  disk  barely  visible,  and  4)  heavy  overcast  solar  disk  not  visible.  The  plotted 
values  are  normalized  and  represent  the  fraction  per  micrometer  of  the  total  incident  ir¬ 
radiance  and  thus  must  be  multiplied  by  the  total  incident  irradiance  from  400-2400  nm.  One 
somewhat  disturbing  feature  of  the  curves  is  that  the  rrt  axima  are  about  700  nm,  which  is  con¬ 
trary  to  other  observations  (Gast  1960)  and  to  our  belief  of  a  yellow  sun.  In  actuality  the  maxi¬ 
ma  are  not  indicative  of  a  reddening  of  the  sun,  but  rather  are  a  manifestation  of  10-15%  un¬ 
certainties  in  absolute  spectral  irradiances.*  This  illustrates  a  need  for  more  precise  meas¬ 
urements  of  incident  spectral  irradiance  under  arctic  conditions. 

Both  the  albedo  and  the  upwelling  irradiance  at  the  surface  can  nave  a  specular  reflection 
component  (R0).  This  is  strictly  a  surface  phenomenon  arising  from  index  of  refraction  dif¬ 
ferences  at  the  interface  between  two  media.  The  indices  of  refraction  of  ice  and  water  are  so 
close,  1.31  and  1.33  respectively,  that  Rlt  at  the  bottom  ice-water  boundary  is  assumed  to  be 
zero.  However  since  air  has  an  index  of  refraction  of  1,  R()  at  the  surface  is  not  negligible. 
Under  diffuse  sky  conditions,  R0  for  air/ ice  or  air/ water  will  be  0.05  (Perovich  and  Grenfell 
1982).  The  dry  and  wet  snow  types  have  a  rough  granular  surface  and  we  cannot  distinguish 
specular  reflection  from  scattering,  so  we  assume  Rn=  0. 


'  Personal  communication  with  T.C.  Grenfell,  University  of  Washington,  1988. 


6 


A  computer  program,  utilizing  thesolution  presented  in  the  theory  section,  'vas  written  to 
calculate  the  upwelling  and  dowmvelling  irradiances  and  associated  par-vi.'ters  for  an  ar¬ 
bitrary  medium  consisting  of  up  to  15  distinct  layers.  By  inputting  the  reflection  and  ex¬ 
tinction  coefficients  for  the  nine  media  types,  the  general  model  can  be  customized  for  sea 
ice.  l'lie  model  allows  selection  of  the  number,  type,  and  thickness  of  layers,  along  with  the 
incident  shortwave  irradiance,  shv  conditions  and  spectral  range  of  interest.  Output 
parameters  include  spectral  values  of  surface  up  welling  irradiance,  bottom  down  welling  ir¬ 
radiance.  albedo,  transmittance,  and  PAR,  and  wavelength  integrated  values  of  bulk  albedo, 
transmittance,  transmitted  flux,  and  PAR.  Because  of  the  relative  simplicity  of  the  model, the 
solution  can  be  quickly  calculated  on  a  personal  computer. 


APPLICATIONS  OF  THE  MODEL 

We  shall  demonstrate  the  utility  of  the  radiative  transfer  model  by  applying  it  to  three 
distinct  cases.  1  >  obtaining  a  two-dimensional  spatial  picture  of  the  transmitted  radiation 
field  through  a  cold  ice  cover  which  varies  in  X  and  V,  2)  monitoring  temporal  changes  in 
light  transmission  through  melting  ice,  and  3)  determining  combined  spatial  and  temporal 
variations  m  light  transmission.  The  major  locus  of  these  three  examples  will  be  on  the  trans¬ 
mitted  spectral  radiation,  with  some  attention  paid  to  spectral  albedos.  In  addition  to  show¬ 
ing  the  capabilities  of  the  model,  the  examples  have  been  selected  to  illustrate  the  effects  of 
\  anations  in  ice  thickness,  snow  cover,  surface  conditions,  and  ice  types  on  the  optical  prop¬ 
erties. 

Spatial  variations 

Tucker*  used  a  hot  water  drill  to  perform  a  detailed  mapping  of  ice  and  snow  thickness  on 
several  floes  in  the  Beaufort  Sea.  Figure  5  summarizes  the  snow  and  ice  thickness  data  from 
one  of  these  floes,  where  121  ice  thicknesses  were  measured  for  a  100-x  100-m  section  using 
a  10- m  grid.  This  area  consisted  of  deformed  multiyear  ice  with  a  primarily  columnar  crystal 
texture.  Since  the  observations  were  made  in  April, the  snow  cover  was  cold  and  dry,  and  the 
ice  was  cold  but  not  below  the  eutectic  point.  Ice  thicknesses  ranged  from  1.8  to  6.5  m  with 
a  mean  of  3.01  in.  while  snow  thicknesses  varied  from  0  to  0.34  m  with  a  mean  of  0.07  m. 

The  ice  cover  is  assumed  to  be  relatively  homogeneous,  consisting  of  interior  white  ice. 
Thus  variations  in  light  transmission  are  primarily  due  to  differences  in  snow  and  ice  thick¬ 
ness  The  two-dimensional  spatial  distribution  of  transmitted  light  at  470  nm  is  plotted  in 
Figure  n.  The  fraction  transmitted  spans  over  5  orders  of  magnitude  from  0.88  x  10""  to  0.63 
.<  10  .  Features  in  this  plot  do  not  necessarily  correlate  with  the  ice  thickness  topography.  The 
peaks  in  transmitted  light  correspond  to  areas  of  bare  ice,  while  minimum  values  occur 
where  the  snow  was  the  deepest.  Distributions  at  other  wavelengths  are  similar  in  shape  to 
Figure  6,  but  have  reduced  magnitudes. 

The  question  arises  whether  this  detailed  spatial  computation  results  in  a  better  represen¬ 
tation  of  light  reflection  and  transmission  for  this  area.  To  address  this  issue,  spectral  albedos 
and  transmittances  were  calculated  for  the  121  points  on  Tucker's  thickness  grid.  These 
values  were  then  averaged  to  generate  an  areal  estimate  for  the  region  and  are  plotted  as 
curve  1  i  n  Figures  7a  a  nd  7b  For  comparison  another  areal  estima  te  was  computed  more  sim- 
plisticallv  bv  first  determining  the  mean  snow  and  ice  thicknesses  and  then  calculating  spec¬ 
tral  values  for  this  single  case  (cur  e  2).  As  curves  1  and  2  indicate,  there  is  a  significant  dis¬ 
parity  between  these  two  met  hex.  .  ith  the  simplistic  method  yielding  albedos  roughly  0.05 
higher  and  transmittances  lie  i-  r  order  of  magnitude  lower  than  the  more  detailed  anal¬ 
ysis.  Of  course  this  is  not  sui  p  .  since  radiative  transfer  processes  are  more  exponential 
than  linear  with  depth.  Asa  di.ect  uence  of  this,  the  actual  areal  averaged  transmittance 

will  always  be  greater  than  or  eq  ,.  to  values  computed  by  using  mean  snow  and  ice 


*  IVison.il  c  omnumuMtion  with  Walter  B.  Tucl  er  III,  CRRELJMSS. 
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Figure  5.  Snow  and  ice  thickness  data  from  Tucker  (personal  communi¬ 
cation, 1988).  The  bottom  and  top  surfaces  represent  the  thicknesses  of  the 
ice  and  the  snow  cover  respectively.  The  thickness  data  were  taken  at  1 0- 
m  spacings  on  a  100-  x  100-m  grid. 


thicknesses  for  the  area.  The  difference  is  particularly  pronounced  for  cases  such  as  this  in 
which  part  of  the  area  is  snow  free. 

Temporal  variations 

Even  a  spatially  uniform  ice  sheet  can  be  optically  complex  when  it  is  undergoing  temporal 
changes.  From  an  optical  perspective  the  most  dramatic  changes  occur  in  the  early  portion 
of  the  melt  season  as  the  ice  cover  evolves  from  an  opaque  snow-covered  medium  to  trans¬ 
lucent  bare  or  ponded  ice.  This  transition  is  of  particular  interest  biologically,  since  some  of 
the  organisms  living  in  and  under  the  ice  are  light  limited  and  the  removal  of  the  snow  cover 
governs  the  onset  of  biological  activity.  To  illustrate  this  point  we  will  examine  the  optical 
evolution  of  a  topographically  simple  ice  sheet  during  the  early  melt  season.  In  the  spring  of 
1982,  Mould  Bay  fiord  (Prince  Patrick  Island  in  the  Canadian  Archipelago)  was  covered  by 
2.2  m  of  flat,  undeformed,  first-year  ice  overlaid  by  0.30  m  of  snow  (Perovich  1983).  The  snow 
melted  quickly  between  21  and  25  June,  transforming  the  entire  ice  sheet  into  a  melt  pond 
for  a  few  days,  until  it  isostatically  readjusted,  leaving  predominantly  melting  blue  ice  and 
shallow  melt  ponds,  with  some  patches  of  drained  white  ice.  These  changes  in  surface  con¬ 
ditions  are  illustrated  in  Figure  8.  By  10  July  the  ice  had  thinned  to  1.46  m,  with  most  of  the 
melting  occurring  during  July.  Table  2  summarizes  ice  thickness,  snow  depth,  pond  depth. 
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Figure  6.  Spatial  distribution  of  light  transmittance  at  470  nm  for  ice 
and  snow  conditions  of  Figure  5. 


Wavelength  (^1  Wavelength  (Mm) 


Figure  7.  Areally  averaged  values  of  spectral  albedo  and  transmittance  for 
ice  cover  of  Figure  5.  Curve  1  is  the  average  of  the  albedos  and  transmittances 
computed  for  the  121  grid  points.  In  curve  2  the  average  snow  and  ice 
thicknesses  of the  grid  were  determined;  then  using  these  average  thicknesses 
the  spectral  albedos  and  transmittances  were  calculated. 

incident  shortwave  irradiance,  and  cloud  cover  for  20  June  to  10  July  1982.  With  these  data 
as  input,  the  radiative  transfer  model  was  run  for  each  day. 

In  Figure  9  albedo  is  plotted  as  a  function  of  wavelength  and  time.  The  first  decrease  in  al¬ 
bedo  occurred  on  21  June  when  the  snow  began  to  melt.  This  drop  was  somewhat  greater  at 
longer  wavelengths.  There  was  a  precipitous  drop  three  days  later  when  the  snow  melted 
and  the  surface  became  completely  covered  by  melt  ponds.  After  this  point  spectral  albedos 
were  fairly  constant,  with  some  small  variations  due  to  changes  in  pond  depths. 


a.  21  liuie,  inciting  snow.  b.  2  July,  heavily  ponded. 

Figure  S.  Photographs  illustrating  surface  conditions  on  Mould  Bay  ice  cover. 


Table  2.  Environmental  conditions  for  Mould  Bay  ice  cover. 


Date 

Ice 

thickness 

Surface 

conditions 

Incident 
short  lease 

Cloud 

(1982) 

(in) 

(m) 

(W/nr) 

conditions 

20  June 

2.22 

0.30  DS 

310 

Heavy  overcast,  solar  disk  not  visible 

21  June 

i  -n 

0.30  MS 

316 

Light  clouds,  solar  disk  clearly  visible- 

22  June 

2.22 

0.20  MS 

357 

clear  skies 

23  June- 

2.22 

0.10  MS 

352 

Clear  skies 

24  June- 

2.22 

0.10  MP 

209 

Heavy  overcast,  solar  disk  not  Visible- 

25  June- 

2.20 

0.20  MP 

199 

Heavy  overcast,  solar  disk  not  visible 

26  June- 

2.18 

0.10  MP 

178 

Heavy  overcast,  solar  elisk  not  visible 

27  June- 

2.16 

0.05  MP 

353 

Light  clouds,  solar  disk  clearly  visible 

28  June- 

2.14 

0.05  MP 

235 

Cloudy,  solar  disk  barely  visible 

29  June 

2.12 

0.05  MP 

332 

Clear  skies 

30  June 

2.10 

0.1  el  MP 

326 

Clear  skies 

1  July 

2.03 

0.05  MP 

231 

Cloudy,  solar  disk  barely  visible- 

2  July 

1.96 

0.10  MP 

337 

dear  skies 

3  July 

1.89 

0.05  MP 

228 

Cloudy,  solar  disk  barely  visible 

4  July 

1.83 

0.05  MP 

126 

Heavy  overcast,  solar  disk  not  visible 

5  July 

1.76 

0.05  MP 

170 

Heavy  overcast,  solar  disk  not  visible 

6  July 

1.69 

0.05  MP 

241 

Cloudy,  solar  disk  barely  visible 

7  July 

1.62 

0.05  MP 

313 

Light  clouds,  solar  disk  clearly  visible 

3  July 

1.56 

0.05  MP 

335 

Clear  skies 

9  July 

1.53 

0.05  MP 

312 

Clear  skies 

10  July 

1.49 

0.05  MP 

312 

Clear  skies 

Note:  The  surface  conditions  give  the  thickness  and  type  of  surface  layer  where  DS  is 
dry  snow,  MS  is  melting  snow,  and  MP  is  melt  pond. 


Figure  9.  Temporal  evolution  of  spectral  albedo  for  the  Mould  Bay  ice  cover. 


Spectral  values  of  transmitted  PAR  for  the  same  time  period  are  plotted  in  Figure  10.  The 
snow  cover  was  still  present  for  the  first  few  days  and  an  insignificant  amount  of  light  was 
transmitted.  Once  the  snow  melted  and  ponds  formed,  PAR  values  increased  sharply, 
especially  at  shorter  wavelengths.  As  the  ice  thinned  over  the  next  two  weeks  there  was  a 
general  upward  trend  in  transmitted  PAR.  This  upward  trend  increased  in  July  as  the  ice 
thinning  accelerated.  The  occasional  small  decrease  occurred  on  cloudy  days  with  low 
incident  solar  irradiance.  The  spectral  peak  of  transmitted  PAR  was  at  530  nm.  Total  PAR 
displayed  a  similar  behavior,  increasing  by  more  than  an  order  of  magnitude  as  the  the  snow 
cover  began  to  melt  (20-21  June),  then  sharply  by  a  factor  of  100  as  the  snow  completely 
melted  and  the  ice  surface  became  flooded  (21—24  June),  and  finally  a  gradual  10-fold 
increase  as  melting  progressed  and  the  ice  cover  thinned  (25  June-10  July). 

Spatial  and  temporal  variations 

In  the  more  dynamically  active  regions  of  the  Arctic  simple  undeformed  ice  sheets,  as 
found  in  Mould  Bay,  are  quite  scarce.  More  common  is  an  ice  cover  with  some  surface  topog¬ 
raphy,  perhaps  a  ridge,  a  few  hummocks,  a  low  area  or  two,  and  typically  a  variable  snow 
cover.  In  such  a  case,  spatial  variations  in  ice  conditions  can  further  complicate  the  temporal 
evolution  of  the  transmitted  radiation  field. 

To  illustrate  the  combined  effects  of  spatial  and  temporal  variations  let  us  examine  three 
sites  that  were  monitored  from  20  June  to  7  July  during  the  1984  Marginal  Ice  Zone  Experi¬ 
ment  (Maykut  and  Perovich  1985).  These  sites  were  all  within  5  m  of  one  another  along  a 
small  pressure  ridge.  Site  1  was  on  the  side  of  the  ridge  where  the  snow  cover  was  deepest 
and  had  snow  present  through  7  July.  Site  2  was  at  the  top  of  the  ridge,  where  the  surface  was 
snow  free,  though  covered  with  0.03  m  of  melting  granular  white  ice  throughout  the  experi¬ 
ment.  Site  3  was  initially  covered  by  0.21  m  of  snow,  which  melted  by  3  July.  A  melt  pond 
formed  at  this  site,  reaching  a  depth  of  0.20  m  by  7  July.  Snow  and  ice  thicknesses,  surface 
conditions  and  the  incident  radiation  field  for  these  three  sites  are  summarized  in  Table  3. 

Spectral  albedos  for  the  three  sites  are  summarized  in  Figure  1 1 .  Surface  conditions  at  site 
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Figure  1 0.  Temporal  evolution  of  spectral  and  total  PAR  for  the  Mould  Bay  ice 
cover.  Units  of  PAR  are  pE/nrs. 

2  (dashed  line)  remained  constant  during  the  experiment  and  this  was  reflected  in  an  albedo 
that  was  temporally  invariant.  Site  1  (shaded  envelope)  exhibited  a  very  slight  decrease  with 
time  due  to  the  thinning  of  the  snow  cover.  Changes  in  spectral  albedo  were  most  dramatic 
at  site  3,  where  surface  conditions  evolved  from  cold  dry  snow,  to  melting  snow,  to  bare 
white  ice,  to  melting  blue  ice,  to  a  melt 
pond  with  a  corresponding  continual  de-  1  ( 
crease  in  albedo.  Daily  albedos  at  site  3 
from  20  June  to  7  July  are  represented  by 
the  family  of  curves  in  Figure  11.  0-1 

Transmitted  values  of  total  PAR  for  the 
three  sites  are  presented  in  Figure  12.  Total  °  O.i 
transmitted  PAR  integrated  over  the  18-  |$j 

day  period  was  17.6  |iE/m2s  for  site  1,12.6  <  0  , 

for  site  2,  and  31 .6  for  site  3.  Values  at  sites 
1  and  3  were  small  initially,  but  increased  , 
as  the  snow  cover  melted.  Transmitted 
PAR  was  the  greatest  at  site  3  after  the 
snow  melted  and  a  melt  pond  formed.  At 
site  2  results  were  fairly  constant  during 
the  entire  observation  period.  There  was  a 
sharp  increase  in  total  transmitted  PAR  at  Figure  2 2 .  Spectral  albedos  from  sites  l,  2,  and  3.  The 
site  1  between  6  and  7  July  as  the  snow  shaded  area  represents  albedos  for  site  2,  the  dashed 
thickness  decreased  from  0.10  to  0.04  m.  line  shows  the  albedo  at  site  2,  and  the  series  of  curves 

On  20  June  site  2  had  the  highest  trans-  depicts  daily  albedos  at  site  3  from  20  June  (high)  to 
mission,  yet  by  the  end  of  the  experiment  7  July  (low) 


Stte  3  June  20 


400  600  800  1000 
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Table  3.  Environmental  conditions  for  MIZEX  84  sites. 


Date 

Cloinl 

Incident 

slwrtivmv 

Site 

Snmv 

1 

1  sJ  and 

ice  ii)  thicki 
Site  2 

*>>  (111) 

Site 

3 

<19841 

uiMifitl.ms 

(  VV/mi;J 

H 

H 

H 

H 

H 

H 

;0 June 

Clear  skies 

300 

0.38 

2.02 

0.00 

3.87 

0.21 

3.32 

21  Juno 

Clear  skies 

303 

0.34 

2.01 

0.00 

3.83 

0.1S 

3.29 

22  lune 

Light  clouds, 

215 

0.32 

2.01 

0.00 

3.82 

0.16 

3.28 

23  June 

solar  disk  visible 

Light  clouds. 

220 

0.29 

2.00 

0.00 

3.79 

0.15 

3.27 

24  June 

solar  disk  visible 
Clear  skies 

277 

0.2S 

1 .98 

0.00 

3.79 

0.14 

3.26 

23  lune 

Clou  Jv,  solar  disk 

183 

0.27 

1.96 

0.00 

3.78 

0.14 

3.25 

2h  June 

ba:elv  visible 

Cloudy,  solar  disk 

194 

0.25 

1.92 

0.00 

3.7t> 

0.11 

3.25 

27  June 

barely  visible 

Cloudy,  solar  disk 

172 

0.24 

1.85 

0.00 

3.74 

0.11 

3 

28  June 

barely  visible 

Cloudy,  solar  disk 

171 

0.24 

J  .S3 

0.00 

3.73 

0.10 

3.20 

29  June 

barely  visible 

Cloudy,  solar  disk 

159 

0.23 

1.79 

0.00 

3.69 

0.06 

3.16 

30  June 

bareiy  visible 

Cloudy,  solar  disk 

168 

0.19 

1.76 

0.00 

3.65 

0.03 

3.13 

1  July 

barely  visible 

Heavy  overcast. 

114 

0.17 

1.71 

0.00 

3.60 

0.02 

3.11 

2  July 

solar  disk  not  visible 
Cloudy,  solar  disk 

170 

0.16 

1.68 

0.00 

3.59 

0.01 

3.10 

3  July 

barely  visible 

Cloudy,  solar  disk 

157 

0.15 

1.67 

0.00 

3.57 

0.00 

3.06* 

4  July 

barely  visible 

Heavy  overcast, 

136 

0.13 

1.65 

0.00 

3.54 

0.00 

3.00 

3  July 

solar  disk  not  visible 
Heavy  overcast, 

137 

0.11 

1.62 

0.00 

3.49 

0.00 

2.91 

6  July 

solar  disk  not  visible 
Heavv  overcast, 

104 

0.10 

1.60 

0.00 

3.40 

0.00 

2.78 

7  July 

solar  disk  not  visible 
Heavv  overcast 

100 

0.04 

1.60 

0.00 

3.23 

0.00 

2.64 

solar  disk  not  visible 

*  H  for  site  3  designates  the  formation  of  a  melt  pond. 
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Figure  12.  Total  PAR  (\iEjurs)  as  a  function  of  time  for  sites  1,2,  and  3. 


Wav«kingth  (nm) 


a.  Site  1 . 


b.  Site  2. 


Figure  13.  Contours  of  spectral  PAR  for  sites  1,  2,  ami  3. 
The  contours  represent  orders  of  magnitude  of  spectral 
transmitted  PAR  (pE//m2s),  for  example  the  -2  contour 
denotes  a  PAR  of  10  2  pE/nrs.  The  contour  interval  is  0.2 . 


it  was  the  lowest.  In  fact,  readings  at  site  3  increased 
from  20%  of  site  2  to  500%.  of  site  2.  Even  over  this 
small  5-m  spatial  scale,  at  a  particular  time,  differences 
in  transmitted  PAR  between  the  three  sites  were  as 
large  as  a  factor  of  five,  illustrating  the  potential 
inaccuracies  of  estimating  transmitted  PAR  from  a 
single  set  of  measurements. 

Transmitted  PAR  as  a  function  of  time  and  wave¬ 
length  for  sites  1,2,  and  3  are  plotted  in  Figure  13.  The 
contours  are  isopleths  of  the  base  10  logarithm  of 
transmitted  PAR  (0.2  intervals).  The  sharp  dropoff  of 
PAR  at  wavelengths  beyond  650  nm  is  evident  in  all 
three  cases.  Maximum  PAR  is  transmitted  between 
450  and  550  nm.  Temporal  changes  were  greatest  at  site  3,  with  an  increase  in  PAR  at 
of  one  and  a  half  orders  of  magnitude.  In  essence  the  formation  of  the  melt  pond  at 
analogous  to  opening  a  window  in  the  ice  cover. 


470  nm 
site  3  is 


DISCUSSION 

A  sea  ice  cover  can  have  tremendous  variability  in  the  reflected  and  transmitted  radiation 
fields.  On  a  single  floe  ice  thickness  can  vary  from  centimeter-thick  nilas  to  ridges  tens  of 
meters  thick,  while  snow  depths  can  range  from  zero  to  a  meter  or  more.  Over  horizontal  dis¬ 
tances  of  a  few  meters  surface  conditions  can  change  from  snow-covered  ice  to  bare  white 
ice  to  ponded  ice.  As  we  have  seen  in  the  three  case  studies,  this  can  result  in  values  of  light 
transmittance  that  vary  by  several  orders  of  magnitude.  Temporal  evolution  during  the  melt 
cycle  is  also  quite  significant.  Removal  of  the  snow  cover  results  in  a  reduction  in  albedo  and 
a  sharp  increase  in  transmittance.  The  formation  of  melt  ponds  causes  a  further  decrease  in 
albedo  and  enhancement  of  transmittance. 
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This  complexity  makes  estimation  of  areally  or  temporally  averaged  values  of  reflected 
or  transmitted  light  difficult.  However,  such  values  are  important  for  a  number  of  applications, 
including  assessments  of  biological  activity,  internal  ice  melting,  and  regional  albedos. 
Single  point  measurements  of  albedo  or  transmittance  are  inadequate  to  define  areal  aver¬ 
ages,  while  measurements  on  a  detailed  spatial  scale  are  not  practical. 

The  model  presented  in  this  report  provides  a  mechanism  for  estimating  light  reflection 
and  transmission  tor  such  complex  ice  covers.  Because  of  the  model's  computational  sim¬ 
plicity,  the  optical  properties  of  an  ice  cover  can  be  represented  in  spatial  or  temporal  detail. 
The  input  parameters  needed  for  the  computation  are  all  easily  observable  quantities:  ice 
thickness,  snow  depth,  ice  surface  conditions,  and  a  qualitative  description  of  the  medium's 
composition  in  terms  of  ice  and  snow  types.  Once  the  structural  description  is  specified,  the 
modei  enables  spectral  albedos,  PAR  values,  and  transinittances  to  be  calculated. 

Currently  the  model  is  limited  to  nine  snow  and  ice  tvpes.  While  these  types  constitute  a 
considerable  collection,  there  are  a  few  significant  omissions.  For  example,  there  is  a  scarcity 
of  data  on  the  optical  properties  of  antarctic  pack  ice.  Because  of  the  significant  biological 
activity  within  and  beneath  the  antarctic  ice  cover,  this  is  a  region  of  extreme  interest  for  op¬ 
tical  modeling.  Structurally  antarctic  sea  ice  is  quite  different  from  arctic  ice  (Gow  et  al.  1982), 
with  large  expanses  of  pancake  ice  and  a  preponderance  of  frazil  rather  than  columnar  crys¬ 
tals.  This  difference  in  ice  structure  presumably  results  in  some  d  isparitv  in  optical  properties. 
I  he  presence  of  contaminants  in  the  ice,  such  as  biological  organisms,  would  also  change  its 
optical  properties.  Measurements  are  needed  to  determine  spectral  scattering  and  extinction 
coefficients  tor  these  ice  tvpes  to  extend  the  model  to  these  cases. 
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APPENDIX  A:  DECOUPLING  OF  FIRST-ORDER  DIFFERENTIAL  EQUATIONS 

The  upwelling  and  dovvmvelling  irradiances  are  defined  by  the  coupled  differential  equa¬ 


tions 

t IF  (z,X)  =  -k.F.(z,X)dz-r.Fl(z,X)dz  +  r.F^z,X)dz  (Al) 

dF.(z,X)  =  kf.(z,X)dz  -  rF,(z,X)dz  +  r.F^z,X)dz  (A2) 

In  order  to  solve  for  Ft(z,X  and  Fr(z,U)  these  equations  are  reformulated  as  uncoupled 
second  order  equations.  Equations  Al  and  A2  can  be  rewritten  as 

(d/dz  +  k.  +  r )  F^A)  -  rFr(z,X)  =  0  (A3) 

r.FdzX)  +  ( d/dz  -  k.  -  r. )  F,(z,X)  =  0  (A4) 

respectively.  Operate  on  eq  A3  by  multiplying  by  (r.): 

r.(d/dz  +  k.  +  r. )Fl(z,  X)  -  r. 2F.(z, X)  =  0.  (A5) 

Operate  on  eq  A4  with  ( d/dz  +  kt  +  r.): 

(d/dz  +  k  +  r.  )r.  F.(z,X)  +  ( d/dz  +  k  +  r  )  [d/dz  -  (k,  +  r )]  F^X)  =  0.  (A6) 

Subtracting  eq  A5  from  A6  gives 

UF/dz^F^XHk)2  +  2k.r.f)F,(z,X)  =  0 ,  (A  7) 

which  has  the  solution 

Fr(z,X)  =  C  sinhf^z)  +  D  coshftq;:).  (A8) 


An  equation  for  FL(z,X)  is  determined  in  a  similar  fashion.  Equation  A3  is  operated  on  with 
(d/dz-k,  -  r) 

and  eq  A4  is  multiplied  by  r.  giving 

(d/dz  -  k.  -  rK)(d/dz  +  k.+r.)  F.(z,X)  -  (d/dz  - 1  -  r. )  r./T(z,X)  =  0  ( A9) 

and 

r.^F^X)  +  (d/dz  -t-  r. )  r.JFr(z,X)  =  0.  ( A10) 

Adding  eq  A9  and  A10  gives 

(dd/dz2)  FXz,X)  -  (k.+r)  (k/2+2k.rJFl(z,X)  =  0,  (Al  1 ) 

which  has  the  solution 

Ft(z,X)  =  A  sinhfic^z)  +  B  cosh(iqz).  (A12) 

Equations  A8  and  A12  can  also  be  expressed  in  terms  of  exp  (to::)  and  exp[-  (kxz)].  A  solution 
using  sinh(K.  z)  and  cosh(K.z)  was  selected  because  it  gave  greater  numerical  stability  when 
inverting  the  matrix  defined  in  Figure  2. 
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A  facsimile  catalog  card  in  Library  of  Congress  MARC  format  is  reproduced 
below. 
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